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ABSTRACT 

During the past decade several observational and theoretical works provided evidences 
of the binary nature of r/ Carinae. Nevertheless, there is still no direct determination 
of the orbital parameters, and the different current models give contradictory results. 
The orbit is, in general, assumed to coincide with the Homunculus equator although 
the observations are not conclusive. Among all systems, rj Carinae has the advantage 
that it is possible to observe both the direct emission of line transitions in the central 
source and its reflection by the Homunculus, which are dependent on the orbital 
inclination. In this work, we studied the orbital phase dependent hydrogen Paschen 
spectra reflected by the SE lobe of the Homunculus to constrain the orbital parameters 
of ^Car and determine its inclination with respect to the Homunculus axis. Assuming 
that the emission excess is originated in the wind-wind shock region we were able to 
model the latitude dependence of the spectral line profiles. For the first time, we were 
able to estimate the orbital inclination of ryCar with respect to the observer and to 
the Homunculus axis. The best fit occurs for an orbital inclination to the line of sight 
of i ~ 60° ± 10°, and i* ~ 35° ± 10° with respect to the Homunculus axis, indicating 
that the angular momenta of the central object and the orbit are not aligned. We were 
also able to fix the phase angle of conjunction as ~ —40°, showing that the periastron 
passage occurs shortly after conjunction. 
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1 INTRODUCTION 

Among the most massive and luminous stars in the Galaxy 
rj Carinae is possibly the most enigmatic and intriguing ob- 
ject. Unique features, e.g. variable emissions and the out- 
bursts that occured during the 1800's, which generated the 
Homunculus - the giant bipolar gas structure surrounding 
the star -, made difficult the understanding of the real na- 
ture of this object. Its periodicity in the observed light curves 
and spectral lines at all wavelengths indicated its binary na- 
ture with P = 5.52 years (Damineli 1996, Corcoran et al. 
2001, Stahl et al. 2005). 

From the X-ray light curve and spectra it is possible to 
infer the existence of a strong wind-wind shock, expected 
for an O-type or Wolf-Rayet (WR) companion star (Pittard 
& Corcoran 2002). However, since there is no direct obser- 
vation of the individual stars of the system it is difficult 
to precisely derive its orbital parameters, e.g. the eccentric- 
ity, inclination and periastron position, and indirect meth- 
ods have to be used. The periodic variations in the X-ray 
light curve, which is originated by the wind- wind shocked 
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gas during periastron passage, has already provided some 
clues on this issue (Corcoran 2005). Similarly, the sharp 
drop in flux has also been observed at optical (Fernandes- 
Lajus et al. 2008), infrared (Whitelock et al. 2004) and radio 
wavelengths (Duncan & White 2003, Abraham et al. 2005a), 
but with different duration. Corcoran et al. (2001) assumed 
that the "eclipse- like" dips are produced by increased ab- 
sorption in the stellar envelope. In this scenario, the or- 
bital plane should closely intercept the line of sight and 
the periastron should be located in opposition regarding 
the observer. Falceta-Gongalves, Jatenco-Pereira & Abra- 
ham (2005) proposed a different scenario, in which the rem- 
nant post-shocked gas could increase the optical depth and 
explain the observed X-ray light curve if the periastron pas- 
sage is close to conjunction. In a further work, it was also 
possible to reproduce the radio light curves under this as- 
sumption (Abraham et al. 2005b). The same process could 
be invoked to explain the decrease in optical and infrared 
fluxes. 

Although X-rays, optical, infrared, and radio light 
curves may provide clues on the orbit of the system, and 
even be used on the accurate estimation of the stellar wind 
properties, it is well known that modelling spectral line pro- 
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files is the best method to determine orbital parameters of 
binary systems, in general, with good precision. 

Damineli (1996), studying the HeIA10830, reported a 
periodicity of 5.5yrs in the spectroscopic events, the epochs 
in which the line disappear. Steiner & Damineli (2004) 
studying the higher excitation energy lines of He II, on the 
other hand, reported a peak in intensity right before perias- 
tron passage. Compared to the X-ray light curve, the simi- 
larity reveals that they may have the same origin. Further 
studies, with higher resolution spectra (Martin et al. 2006), 
have shown that, not only the line intensities vary within the 
5.5yrs period, but also their profiles. They concluded that 
the line source might be more complex than simply atmo- 
spheric emission, and the decline in flux cannot be explained 
only by eclipsing models. Nielsen et al. (2007) attempted to 
explain the line-profiles assuming that the emission is the 
combination of a major component originated in the stellar 
wind, and other minor contributions at different velocities, 
which origin is not clear. From the line velocity curve, they 
inferred an eccentricity e ~ 0.9, and that periastron should 
be placed near opposition. 

Falceta-Gongalves, Jatenco-Pereira & Abraham (2007) 
showed that line profile variability is expected from emis- 
sion originated in the wind-wind interaction zone. Its study 
showed to be a powerful tool in determining the or- 
bital parameters of WR30a (Falceta-Gongalves, Abraham 
& Jatenco-Pereira 2008), where by fitting synthetic profiles 
to the observational data, it was possible to determine the 
inclination and eccentricity of the orbit. Also, using the line 
velocities obtained directly from the primary O-star, it was 
possible to constrain the stellar masses. 

In the rj Carinae binary system, observations of direct 
and Homunculus reflected emission show multi-peaked spec- 
tral lines, with profiles and intensities that vary with or- 
bital phase. Stahl et al. (2005) showed that the Hell A4686A 
emission lines measured towards the central object and the 
SE lobe of the Homunculus presented different profiles, in 
contradiction to what would be expected assuming spheri- 
cal symmetry. Abraham & Falceta-Gongalves (2007) showed 
that the conical wind-wind shock region could be the main 
source of the observed excess in the Hell A4686A line pro- 
files. To reproduce both the direct and Homunculus reflected 
line profiles it was necessary to assume that the Homuncu- 
lus axis does not coincide with the axis of the orbital plane. 
However, due to the limited amount of data, which were re- 
stricted to orbital phases very close to periastron passage, it 
was not possible to lift the degeneracy regarding the incli- 
nation of the orbital and Homunculus axis. In all cases, the 
observed velocity curve, as well as the line profiles, are well 
reproduced only if periastron occurs near conjunction. 

In the present work we extended the analysis of Abra- 
ham & Falceta-Gongalves (2007) , applying this model to the 
H Paschen 8 line reflected by the Homunculus SE lobe ob- 
served by Stahl et al. (2005) and derive, for the first time, 
the orbital and Homunculus axis inclination. In Section 2 
we describe the model and the adopted geometry for the 
system; in Section 3 we show the results, followed by the 
conclusions in Section 4. 



Orbit axis 



Observer 




Homunculus 
axis 



Periastron 



Figure 1. Scheme of the system geometry. The orbital axis has 
an inclination i regarding the line of sight and i* regarding the 
Homunculus axis. It is assumed that the angle between the line of 
sight and Homunculus axis is of 45° (Davidson et al. 2001, Smith 
et al. 2004). 



2 MODEL 

As the massive and supersonic winds of the two stars collide, 
they create a wind-wind interaction region, structured as 
two shocks at both sides of the contact discontinuity where 
the wind momenta are equal. As described in Luo, McCray 
& Mac- Low (1990), the geometry of the contact discontinu- 
ity is asymptotically conical, with an opening angle /3 given 
by: 



/3 ~ 120° 1 



n 



2/r 



1/3 



(1) 



where rj — M s v s /M p v p , M p and M s are the mass loss rates 
and v p and v s the wind velocities of the primary and sec- 
ondary stars, respectively. 

At the shock surfaces the compressed gas reaches tem- 
peratures of ~ 10 6 — 10 K and emits free-free radiation, 
mostly at X-ray wavelenghts. As this gas flows along the 
shock structure, it cools down to recombination tempera- 
tures and becomes denser (Abraham et al. 2005a,b; Falceta- 
Gongalves, Jatenco-Pereira & Abraham 2005). At this stage, 
it could possibly be the source of He and H lines, explain- 
ing some of the features observed in the spectra of massive 
binary systems (Falceta-Gongalves, Abraham & Jatenco- 
Pereira 2006, 2007). 

The model, proposed in Falceta-Gongalves et al. (2006) , 
considers the emission from each subvolume of the conic 
shaped shock structure between the stars. The line intensity 
is obtained by the integration of each fluid element emission 
along the line of sight: 

^2 " 



I(v) = C(<p) 
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where v a ^ B = ua ow ( — cos /3 cos ip sin i + sin /3 cos a sin ip sin i — 
sin f3 sin a cos i) , is the observed stream velocity component 
projected into the line of sight, a the azimuthal angle of the 
shock cone, i the orbital inclination, ip the orbital phase an- 
gle, a the turbulent velocity amplitude, r the optical depth 
along the line of sight and C(<p) a normalization constant. 

From Eq. 2, it is possible to obtain synthetic line pro- 
files and compare them to the observations, as it was done 
with the Hell A4686A emission line (Abraham & Falceta- 
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Table 1. Parameters and input values 



parameter input value reference 



Gongalves 2007) . To fit the lines reflected by the Homuncu- 
lus, it is also necessary to take into account the angle be- 
tween the Homunculus axis and the line of sight and 
the inclination of the Homunculus axis with respect to the 
orbital plane. 

In Fig. 1 we describe the geometry of the system as we 
assumed in the model. Here, i and i* represent the inclina- 
tion of the LOS and of the Homunculus axis with respect to 
the orbital axis, respectively. In this case, we may write the 
relationship between these angles as: 

cos $ = sin i sin i* cos A0* + cos i cos i* , (3) 

where A<f>* represents the angle between the projection of 
the Homunculus axis into the orbital plane and the line 
of sight. Following Davidson et al.(2001) and Smith et al. 
(2004), we will consider the angle between the Homunculus 
axis and the LOS close to ~ 45°. 



3 OBSERVATIONAL DATA 

Stahl et al. (2005) presented VLT/UVES observations of the 
line emission reflected by the SE lobe of the Homunculus, 
carried out from December 2002 to March 2004, or phase in- 
terval between 0.9 and 1.127, including the minimum event 
of 2003.5. The Paschen P8 line profile presented strong vari- 
ability, with clear P Cygni absorption, probably related to 
the stellar wind. Also, in all observations, two prominent 
excess emission bumps are visible, which we interpreted as 
extra emission from the shock region. The digitalized excess 
emission velocity profiles are shown in Figure 2 as filled cir- 
cles, corrected by the expansion velocity of the Homunculus 
V = 100 km s _1 . The numbers at each spectrum represent 
the date of the observation (JD - 2.400.000). 

The observed peaks reflected by the Homunculus pole 
occurred at v ~ 0km s -1 and v ~ +220km s _1 , while direct 
observations towards the central star revealed v ~ and v < 
150km s _1 (Weis et al. 2005). This difference shows that, at 
least part of the observed emission has an anisotropic origin. 
We assumed that this latitude dependence is connected to 
the emission from the conic shock region. 

4 RESULTS 

We applied the model described in the Section 2 to the ob- 
servational data. For the wind parameters we used M p = 



2.5 x 10~ 3 M yr" 1 , v p = 700 km s" 1 , M s = 5 x 10~ 5 M 
yr _1 , v s = 3000 km s _1 , which results in r\ = 0.1 and, from 
Eq. 1, we obtain (5 — 45°. For the flow along the conic sur- 
face we assumed va ow sin (i = 450 km sT 1 and a — 0.33wfl ow 
(Abraham & Falceta-Gongalves 2007) . Because of the size of 
the Homunculus, we must also take into account the delay 
in the light travel from the wind-wind interaction zone to 
the reflection nebula. The data obtained from the reflection 
of the expanding nebula corresponds to an emission that 
occurred in a phase previous than that observed directly 
towards the central star. This delay is particularly impor- 
tant for this model because a difference in the phase of the 
emission corresponds in a difference in the orientation of 
the shock cone regarding the Homunculus and, therefore, a 
difference shift velocity. The time delay may be calculated 
as: 

^dclay ^expansion 
C ^expansion 

considering an expansion of 100 km s _1 since the outburst 
about 160 years ago, we get ideiay ~ 20 days. Therefore, in 
order to precisely fit the model, we must subtract 20 days in 
the orbital phase and determine the actual orientation of the 
shock cone for the given phase angle. For the determination 
of the orbital phases we used the epoch of conjunction as 
June 29, 2003 (Abraham et al. 2005b), e = 0.95 and an angle 
between periastron and conjunction 0o = 40° (Abraham & 
Falceta-Gongalves 2007) 

In Table 1 we show the input parameters of the model 
and their corresponding values. In the present work the in- 
put parameters, shown in Table 1, are fixed. We will briefly 
discuss the role of each of them on the obtained results later 
in the paper. Finally, we fitted the observational data us- 
ing Eq. 2, together with Eq. 3, for an optically thin gas 
(r ~ 0). In Falceta-Gongalves et al. (2007) it has been 
shown that optically thick wind-wind shocks result in sin- 
gle peaked profiles, while the data of Paschen lines of r\ Car 
presented double-peaked profiles at all phases. In order to 
obtain the best fitting set of values of i and i*, we calculate 
the synthetic line profiles for the observed epochs for the 
range —180° to 180°, independently for i and i* , in steps of 
Ai — 10°. The observed data and three of the best mod- 
els we fitted are shown in Fig. 2. The lines represent the 
models for different values of i and i* . It is noticeable from 
the modeled spectra that the line profiles are very sensi- 
tive to changes in the orbital inclination regarding both the 
observer and the Homunculus. The best fit, i.e. the mini- 
mum cumulative error for all orbital phases, is obtained for 
i* ~ 35° ±10° and i ~ 60° ±10°, i.e. an angle between the or- 
bital plane and the line of sight of ~ 30° , in agreement with 
the radio detection of the almost edge-on gas disk formed 
around the system. Note that the error of Ai = 10° repre- 
sents the "resolution" of the sets of i and i* that we have 
calculated. Due to a small degeneracy between i and i* for 
small values of Ai, it is not possible to reduce the impreci- 
sion in the determination of these parameters. 

The most important issue regarding these results is that 
the Homunculus and the orbital plane axis are not aligned, 
as usually assumed (see Steiner & Damineli 2004, Martin 
et al. 2006, Okazaki et al 2008). Since the Homunculus is 
assumed to be coincident with the stellar rotation axis at 
the epoch when the outburst event occurred, we may con- 
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Figure 2. Time series of the Paschen 8 line profiles, after continuum subtraction, measured for the Homunculus SE lobe. Left: the lines 
show the modeled line profiles for i* = 35° and orbital inclinations regarding the line of sight i = 30° (dashed), 60° (solid) and 80° 
(dotted). Right: the lines show the modeled line profiles for i = 60° and orbital inclinations regarding the Homunculus axis i* = 55° 
(dashed), 35° (solid) and 15° (dotted). 



elude that the stellar and orbital momenta were not parallel 
(or even close to that). Considering other massive binary 
systems, it is well known that the relaxation timescales for 
close binary systems are short, of order of hundred to a few 
thousands of years (e.g. Tassoul 1990). Should we expect the 
orbital and stellar angular momenta to coincide after the 
stellar evolution (~ 10 6 yr) due to tidal relaxation? 



4.1 Relaxation timescales 

Close massive binary systems are known to present strong 
tidal torques (Lecar, Wheeler & McKee 1976) , which result 
in i-) orbital circularization, ii-) stellar rotation synchroniza- 
tion with the orbital period and in-) alignment of stellar and 
orbital angular momenta. The last is particularly important 
considering the results obtained in this work. 

To verify the plausibility of a non-aligned orbital and 
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stellar momenta scenario, we may estimate the relaxing 
timescale of the rj Carinae system. Several authors have ad- 
dressed the problem is the past (e.g. Zahn 1977, Hut 1982, 
Eggleton, Kiseleva & Hut 1998). The gravitational pull force 
acting on each star results in the non-spherical distribution 
of the stellar mass. However, since the stars are moving fast 
in close binaries, the orientation of the tides are not coin- 
cident with the direction of the centers of mass. The ex- 
tra force component is responsible for the net torque that 
changes the energy and angular momentum of the orbital 
motion and stellar rotation. 

In particular, Hut (1982) studied the dynamical evolu- 
tion of close binary systems with high eccentricity due to 
tidal friction between the two stars - a scenario we believe 
is correct for r\ Carinae. In the case of high eccentricity, and 
a ratio of rotation to orbital angular momenta n « 1, the 
orbital and stellar angular momenta tend to become parallel 
for timescales of order: 

™«^o.o2T(°)*(ay i i I ^ (2 _ ()W , (5) 

where k is the apsidal motion constant, D is the distance 
between the stars at periastron passage, q = M s /M p is 
the ratio between the stellar masses, e = (1 — e), r g = 
[Jp/(M p i?p)] 1/ ' 2 , I p is the moment of inertia, S7 P the angular 
velocity and R p the radius of the primary star, 

ojp = [G(Mp + M s )(l + e)] 1/2 D- 3/2 , (6) 

is the orbital angular velocity at periastron, and: 



where r is the constant time lag of the tidal axis with respect 
to a coordinate system corotating with the primary star. The 
apsidal motion constant k depends on the internal structure 
of the star, and is typically of the order of 10 -3 — 1CP 1 . 

From Eqs. 5-7 it is clear that the relaxing timescale 
of the system depends on several poorly known parameters, 
though an estimation is still possible. For D ~ 10R P , r g /q ~ 
1, M s /Mp ~ 0.4, Q p ~ Up and T ~ 10 8 s we have T a u g n ~ 
1Q i6 _ 10 i8 gj or ^ 1Q 8 _ io 10 yr , for k = 1CT 3 - KT 1 . From 
the parameters used in this estimation, the ratio Q. p /u) p is 
the main unknown. Typically, we would expect a ratio ~ 1 
for a relaxed system, but a ratio > 1 for a young wide binary 
system. In this case, the relaxation timescale would be even 
larger and our conclusions unchanged. 

Typically, the lifetime of stars with M > 60M Q is 1 — 5 
Myr (Schaller et al. 1992). Therefore, we may conclude that 
7) Carinae system, during its evolution up to the outburst 
in the 19th century, had not enough time to relax, and we 
should not expect the orbital and rotation angular momenta 
to be aligned. Interestingly, this also reveals that the orbital 
high eccentricity (e > 0.9) may be primordial, even though 
it could also be a result of the massive mass ejection events. 



5 CONCLUSIONS 

Although many authors agree about the binary nature of 
t] Carinae, its orbital parameters and orientation are still 
poorly known. From the X-ray and radio light curves, as well 
as the spectral line velocity curves, it is possible to derive 



high eccentricity (e > 0.9) (Ishibashi et al. 1999, Abraham 
et al. 2005b, Falceta-Gongalves, Jatenco-Pereira & Abraham 
2005). Also, from the spectral line profiles it was believed 
that the periastron should happen in opposition to the ob- 
server, with 0o = 270° , if one assumes that the emission 
comes exclusively from the stellar winds (Davidson 1997, 
Nielsen et al. 2007). Abraham & Falceta-Gongalves (2007) 
showed that the wind-wind shock region can be responsible 
for the HeII4686 emission excess, and consistently fit the ob- 
served spectra. However, in this case the periastron should 
occur near conjunction instead (4>o = 40°). 

Regarding the orbital inclination, it is commonly as- 
sumed that the orbital plane should be nearly aligned with 
the Homunculus equator (Davidson 1997, Smith et al. 2004), 
which is inclined ~ 45° with respect to the observer. In con- 
tradiction to that, high resolution radio observations show 
the presence of a gas torus surrounding the system, which 
seems to be seen edge-on (Duncan & White 2003), indicat- 
ing that the orbit may not be aligned with the Homuncu- 
lus. In this scenario, the torus would also be responsible for 
the shell-like effect seen in the X-ray light curve (Falceta- 
Gongalves, Jatenco-Pereira & Abraham 2005). 

In this work we studied the possibility of the orbital 
plane and the Homunculus not being aligned. For that, we 
used the hydrogen Paschen 8 line spectral profiles obtained 
from the SE lobe of the Homunculus. This line is supposed 
to originate in the central system and reflected by the ex- 
panding gas. The observed excess flux presented two peaked 
line profiles, which velocities change with the orbital phases. 
This behaviour is common in massive binary systems where 
the line emission originates in the wind-wind interaction 
zone. The conical shock surface would be resposible for an 
anomalous two peaked profile with velocities changing dur- 
ing the orbital period as the line of sight is intercepted by 
different vector projections of the flowing gas, as detailed in 
Falceta-Gongalves, Abraham & Jatenco-Pereira (2006). 

We applied the synthetic line profile model, consider- 
ing the emission from the central system, and derived the 
reflected fluxes at the SE lobe of the Homunculus. We fixed 
the stellar and orbital parameters f3 = 45°, e = 0.95, 
WflowSin/3 = 450 km s _1 , 0o = —40° and a — 0.33ufl ow 
(Abraham &: Falceta-Gongalves 2007), and left the orbital 
inclination regarding the line of sight (i) and regarding the 
Homunculus axis (i*) as fit parameters. We showed that the 
line profiles are very sensitive to changes in the orbital in- 
clination regarding both the observer and the Homunculus 
axis. Therefore, this method has revealed helpful in con- 
straining the orbital orientation of the system. The best fit 
was obtained for i ~ 60° ± 10° and i* ~ 35° ± 10°. This 
result shows that the line of sight must lie close to the or- 
bital plane (i ~ 30°), as indicated by radio observations and 
theoretical models regarding the X-ray light curve. 

In the presented calculations, most of the orbital and 
stellar parameters were fixed, based on previous references 
(Table 1). It is important though to understand if changes 
in any of these would reflect in major changes on our re- 
sults and conclusions. The eccentricity of the orbit is a well- 
known parameter as it is constrained by the velocity curves 
of several spectral lines - we should expect its value to vary 
within the range e = 0.9 — 0.95. A different eccentricity 
would simply change the orbital phase, in which a given 
synthetic line would fit the data, therefore being irrelevant 
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for the synthetic profiles. The same occurs for <f>o, P, a and 
to- Actually, 0o is also related to the curvature of the shock 
cone. The turbulence amplitude a, on the other hand, is 
responsible for the broadness of the synthetic line profiles. 
Since the observed broadness for each peak is unchanged 
at different orbital phases, its value is also irrelevant on 
the determination of i and i* . It is important though to 
remember that a presents a degeneracy with «fl ow , i.e. if 
the system has a smaller value of «fl ow , we would need a 
larger value of a to account for the same line broadness be- 
ing iifl ow directly obtained from the observed line velocities, 
but it has no influence on the determination of the inclina- 
tion of the orbital plane. Among the fixed parameters, j3 is 
the only strongly correlated with i and i* . As fully discussed 
in Falceta-Gongalves et al. (2006), a sequence of synthetic 
profiles may be well reproduced if an increased value of f3 
is balanced by a decrease of i, regarding the line of sight. 
However, considering the observed profiles, this degeneracy 
is limited to the range in which i < (3 because, otherwise, 
models with i > j3 result in single peaked profiles in certain 
orbital phases. Since the observational data presents two- 
peaked profiles in all phases we may constrain i < (3. Fortu- 
nately, P depends on the stellar winds of the two stars, which 
are approximately well constrained to provide f3 ~ 40° —50°. 
This range of possible f3 values is coherent with the obtained 
range of i ~ 30° ± 10°, regarding the line of sight. 

Also, the results show that the Homunculus axis is not 
perpendicular to the orbital plane. This suggests that the 
stellar angular momentum, probably coincident with the Ho- 
munculus axis, is not aligned with the orbital motion. We 
explained this fact showing that the relaxing timescale for 
momenta alignment and orbital circularization of r\ Carinac 
is longer than the stellar lifetime. This fact also explains the 
high eccentricity of the system (e ~ 0.95), since tidal inter- 
actions did not have enough time to circularize the orbit. It 
is true though that the outburst occured in the 19th century 
could have been responsible for large changes in the orbital 
parameters as well. 

Finally, direct observations toward the central star re- 
vealed an increase in P Cygni absorption during the mini- 
mum (Davidson et al. 2005), as opposed to the polar data, 
which shows almost constant absorption. This fact tell us 
that the main responsible for this absorption is not the pri- 
mary star wind but, if the orbital plane is close to the LOS, 
the secondary star wind or the remaining post-shocked gas 
within the star and the observer. 
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